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ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ 


УДК 539.3:534.1 
А.Н. СОЛОВЬЕВ, А.О. ВАТУЛЬЯН, А.С. СПОЖАКИН, С.Н. ШЕВЦОВ 


РЕКОНСТРУКЦИЯ ДЕФЕКТОВ В СЛОИСТЫХ КОМПОЗИТАХ 


Предлагаются методы реконструкции двух типов дефектов в слоистых композитах: 
расслоений (интерфейсные трещины) и разрывов слоев (поперечные трещины). 
Обратные геометрические задачи реконструкции трещин сводятся к решению 
последовательности систем граничных интегральных уравнений (ГИУ), полученных 
на основе принципа взаимности работ при установившихся колебаниях 
исследуемых образцов. Решение ГИУ проводится на основе сочетания метода 
конечных элементов (МКЭ) метода граничных элементов (МГЭ) и метода 
регуляризации А.Н. Тихонова. Дополнительной информацией для решения 
обратных задач служит поле смещений, измеренное на части внешней границы 
тела, свободной от механических напряжений. Рассмотрены численные примеры 
идентификации расслоений и разрывов в рамках плоской задачи теории упругости. 
Ключевые слова: неразрушающий контроль, реконструкция трещин, композиты, 
граничные интегральные уравнения, метод граничных элементов. 


Введение. Одной из математических моделей деформированного твердого 
тела с дефектами в виде трещин при пренебрежении взаимодействием 
внутренних поверхностей является линейно-упругое тело с разрезами. На 
берегах этих разрезов задаются граничные условия в напряжениях 
(в частности, условия их отсутствия). В рамках такой линейной модели 
правомочно ставить задачи об установившихся колебаниях тела. 
Измерение амплитуд граничных волновых полей позволяет гораздо 
эффективнее, чем в статике, проводить восстановление напряженно- 
деформированного состояния (НДС) внутри тела и на его границах (в том 
числе внутренних), а по структуре этих полей идентифицировать дефекты. 
Среди задач реконструкции трещин внутри упругого тела наиболее 
простыми представляются задачи, в которых известно сечение тела 
(в общем случае криволинейная поверхность), содержащее дефекты. В 
случае плоского сечения его определение для уравнения Лапласа 
проводится в [1], для гармонических колебаний изотропного упругого тела 
в [2]. Методы, разработанные в этих работах, опираются на возможность 
измерить на всей границе тела, как вектор напряжений, так и вектор 
смещений. К задачам такого же типа относятся задачи определения 
интерфейсных трещин на внутренних границах составного упругого тела. 
Некоторые методы решения таких обратных задач теории упругости 
представлены в литературе. Так, в работе [3] предлагается метод 
неклассических ГИУ [4], в [5] применяется итерационный метод, 
основанный на алгоритме, предложенном в [6]. В работе [1] для уравнения 
Лапласа разработана полуявная схема реконструкции системы трещин. 
Одним из существенных требований к постановке этих обратных задач, 
приближающих их к практическому применению, является условие, при 
котором возможно измерение граничных сопряженных полей не на всей 
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внешней границе, а лишь на ее части. В такой постановке в настоящей 
работе предлагается два подхода к решению проблемы идентификации 
расслоений и формулировка систем ГИУ либо относительно скачков 
смещений на трещинах, либо компонент вектора напряжений на 
внутренней границе, содержащей дефекты. Эти системы являются 
интегральными уравнениями Фредгольма 1-го рода с гладкими ядрами, 
поэтому процедура их решения требует регуляризации. В работе решение 
их строится на основе сочетания идей метода граничных элементов, 
метода конечных элементов и метода регуляризации А.Н. Тихонова [7]. 
Полученные результаты опираются на подходы, развитые в работе [8]. 

Постановка обратной задачи. В декартовой прямоугольной системе 


координат Охх,х, (х= (х,х,,х,)) рассматривается конечное составное 
упругое тело, занимающее область Т-= ый (рис.1, слева) и 
ограниченную поверхностью 5. Подобласти Г, представляют собой слои, 
5: - внутренние интерфейсные поверхности раздела подобластей. На 
поверхности 5 задано разбиение 5= 5, 15, 05, (части которого, не 
пересекаются), часть 5, закреплена, на 5, задан вектор напряжений, 5» 
свободна от напряжений. На поверхности 5, выделяется, часть 5, 
свободная от напряжений и доступная для измерения вектора смещений ( 


М 
5, 65,). Имеется система непересекающихся трещин г= Л.Г А 


Г, = Р- в которые расположены на внутренних поверхностях 5, 


(расслоения), рис.1. Краевая задача, в которой кроме определения 
характеристик НДС, требуется определить геометрию системы трещин Г, 
формулируется следующим образом. Она состоит из уравнений линейной 
теории упругости [9] в случае установившихся колебаний: 


(© 2 (®) Е 
б-р, Как, (1) 
К и ®. 
0 й Ст Ив Г (2) 
граничных условий прямой задачи 
| _ ФГ В о. 
и, [5 = 0, 1 $=9; = р), й 5,20; (3) 
условий непрерывности на 5. \Г: 
| | | | 
|) || || || 
ВНИИ (9 
| | | | | | | | | 
| 























условий на берегах трещин: 
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[=0, 4=12,...М, (5) 


и дополнительных условий, отвечающих измерению вектора смещений 
на Э,: 
10 

|5, = и; , (6) 
(®) 


К К 
где 0‘ - Си - Компоненты тензоров напряжении и упругих постоянных; 


“ 


1 


и 


(К 


и - компоненты вектора смещений; р“, @ - плотность и круговая 


частота колебаний; п, - компоненты единичных векторов внешней 
нормали к соответствующим поверхностям. 


0) 
в 5, _ 





Рис. 1. Образец с дефектами и без дефектов 


Замечание 1. На практике обычно задается не распределение вектора 
перемещений всюду на 5, , а значения перемещений в некотором наборе 
точек, соответствующих местам установки датчиков, причем аналог (5) 
имеет вид: 

и (х,)= ии, Х„ 9. (7) 
Вспомогательные задачи 
Задача [. Рассмотрим краевую задачу Т для тела /' без дефектов, тогда на 
внутренней поверхности 5„ граничные условия будут соответствовать 


непрерывности векторов смещений и напряжений. 


Эта задача состоит из дифференциальных уравнений движения (1),(2) 
относительно и“”, х И, ‚ граничных условий 


&)* &)* й&)* 
и |= 0, "= с У, [= 0, (8) 


К ая, = 6.055), 5 5 (9) 


и условий непрерывности на ›.„ : 
(К)* = (К)* — +(+1* 
и; 2 И, Г, | = Г |5 


т Е" 








ри 5? (10) 


Задача П. Далее рассмотрим краевую задачу П для тела 
РЕ К < К при действии тех же нагрузок на 5.,, причем 


поверхность в для последнего слоя свободна от напряжений 


ше 


(рассмотрим вначале случай, когда внутри тела У? дефектов нет). 
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Рис. 2. Последовательность вспомогательных задач 
( слева А\ = 1; справа А! = 2) 


Эта задача описывается дифференциальными уравнениями (1),(2) 
относительно й,, х У’, граничными условиями вида 


И, |5. = 0, р |5, = буп, |5 =0, (11) 
ео `В, (12) 

а также условиями отсутствия нагрузок на 5. 
й |120. (13) 


Замечание 2. Зависимость функций О’ (х,6), О, (хЁ) и от 6 означает, 


что может быть рассмотрено некоторое однопараметрическое семейство, 
например 

ОЕ 98. (14) 
Замечание 3. Методология построения ГИУ с помощью задачи П и не 
изменится, если тело у’? будет содержать известные дефекты (возможно, 
найденные на предыдущем шаге). 
Вывод ГИУ с помощью решения задачи Т. Рассмотрим случай наличия 


расслоений и отсутствия разрывов. Предположим, что найдено решение 
классической задачи 1 и обозначим 


в" [5 = 6 °(х,ё ), х В Г 5% * 
Введем в рассмотрение оператор 
(и, р,9`,0°)= шо), - бобр, 45, = 6.) (15) 


5о 5, 


Применив к телам Г и У' теорему взаимности работ [9] и 
учитывая условия непрерывности (4) на 5.„\Г и (10) на 5.„, получим: 


и”. 6945, = 6), ` 


Г 
где 1.0) - скачки компонент вектора перемещений на трещинах, при этом 
на 5.‚„ \Г этот скачки равны нулю. Обозначим: 
4: при хГ 


Х. = 
(®) О при хГ5\ 
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Тогда соотношение (16) представляет собой систему интегральных 
уравнений (при интегрировании по известной поверхности 5) 


относительно функций Х ;(х), х Эш: 


в" (х,6 )Х ‚00а, = <), 6 5%. (17) 


Уравнение (17) при различном выборе частот и фиктивных 
нагрузок может служить для реконструкции трещин, которые являются 


носителями для функций Х ‚ (х). 


Вывод ГИУ с помощью решения задачи 11. Так же, как и в задаче Т, 
рассмотрим случай наличия расслоений и отсутствия разрывов. Пусть 
вспомогательная задача 1 решена, обозначим 


п уе 5. 
Введем в рассмотрение оператор: 
Е(и’,р,д ‚О )= и’ (00, (х6)а5,- 9, (55 )р. 0045, = ЕЕ). (18) 
50 5, г 
Рассмотрим тело, занимающее подобласть У”, и применим к нему 
теорему взаимности работ [9]. Тогда получим, что 
(®) (Ку) хе 
Е и; (с 1, (х)а5, ме ЕС). (19) 
Учитывая, что сомножители подынтегрального выражения в (19) 
определены всюду на 5^' и в силу (5) 


(КК | < 
1, ' | = 0, (20) 
то соотношение (19) может быть представлено как система ГИУ по 
К К. 
известной границе 5: 


(0 (К) & 
Г и, (с 1, (х)45, = ЕС), Я 5. (21) 
При численной реализации решения системы ГИУ (21) в отличие от 
системы (17) основой идентификации трещин служат два свойства 
решения: соотношения (20) и сингулярное поведение компонент вектора 
напряжений на краях трещин. 
Частотное сканирование и регистрация трещин. Если частота 
колебаний рассматриваемой конструкции не определяется 
технологическими условиями, то возможно расширение области 
определения функций, стоящих в левых и правых частях соотношений (17) 
и (21), на следующую частотную область: 


№ 
00 = бе, 


ЕВ п п ' 


где набор интервалов 0 п выбирается из предварительного модального 


анализа конструкции без дефектов, и включает в себя такие формы 
колебаний, при которых происходит интенсивное раскрытие трещин (часто 


это моды растяжения-сжатия или сдвига в окрестности поверхности лу ‚ В 


то время, как изгибные моды менее чувствительны к наличию 
трещиноподобных дефектов). 
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Теперь система ГИУ (17) примет вид: 
(®* с 
‚йе 0 88, = 660), $ 5, 00, (22) 


поскольку носитель функций Х,(х,0) не зависит от @ , то решение 
уравнения (22) в наборе частот из 0 значительно повышает точность 
реконструкции трещин. 

Отметим, что функция С! С ‚0 ) равна тождественно нулю, если в 
теле отсутствуют трещины. Это обстоятельство дает простой способ 
регистрации наличия трещин в теле, причем и в том случае, когда 
поверхность 5, заранее неизвестна. Для этого достаточно выбрать одну 


из схем нагружения в задаче 1, например (9) при © = А и построить 


зависимость: 
6,@)=С:(@,,0), 00 ' 
0 


по виду которой при знании точности измерений поля смещений и; из (6) 


и точности вычисления интегралов (15) можно судить о наличии трещин. 
Аналогом системы (21) служит следующая система: 
(к) (Ку) 
| ‚0 Е“ (х,0 45. = 
В (х,6 1: ' (х,@ )45, Е(@,0), Г 55, 00 , (23) 


где так же, как и в (22), принцип восстановления геометрии Г не зависит 
от частоты @. 


Замечание 4. Отметим, что при выборе функций О, О в рамках 


дискретной схемы их носители выбираются в точках измерения смещений 
х„„ ‚ благодаря чему интегралы в операторах (15), (18) вычисляются явно. 
Численные примеры реконструкции трещин. В качестве первого 
примера применения предложенного метода рассмотрим схему [8], 
приводящую к уравнению (21) в задаче об идентификации двух трещин Г, 
=КЕ и Г,=ММ, находящихся между слоями полупассивного биморфа, 
представляющего собой составной прямоугольник АОВСЕВ (А(0;-0,03), 
0(0;0), В(0;0,03), С(0,1;0,03), Е(0,1;0), 0(0,1;-0,03), К(0,03;0), К0,05,0), 
М(0,065;0), №(0,08;0) - размеры в метрах). Верхний слой ОВСЕ выполнен из 
меди, нижний АОЕР - из пьезокерамики Р7Т-4. Задача решается в условиях 
плоской деформации, колебания (с частотой Л=@ /(2м)=20 кГц) 
возбуждаются разностью потенциалов Г, = 10008, подаваемой на 
электроды, находящиеся на сторонах ОЕ и АБ, сторона АВ жестко 
защемлена, остальные внешние границы свободны от механических 
напряжений. 

При проведении численных экспериментов решение прямой задачи 
и построение ядра уравнения (21) модальный анализ конструкций, 
решение задач для тел без дефектов и "измерение" полей смещений 
производилось с помощью конечно-элементного комплекса АСЕЁГАМ [10]. 
Трещины моделировались отверстиями, у которых поперечный размер 
много меньше их длины (в расчетах их отношение имело порядок 10`”), 
при этом предполагалось, что берега трещины не взаимодействуют между 
собой. На краю трещин конечно-элементная сетка сгущалась с помощью 
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введения дополнительных узлов. В последней версии комплекса АСЕТАМ 
разработан язык команд [11], позволяющий эффективно решать 
множественные задачи, такие как численное построение тензора Грина. 
Считалось, что сторона ВС доступна для измерения вектора смещений и на 
ней было выбрано 39 внутренних равноотстоящих узлов, смещения 
которых моделировали процесс измерений. Частота, на которой 
проводился численный эксперимент, выбрана из соображений 
интенсивного раскрытия трещин на ней при отличии от собственной 
частоты краевой задачи П. В табл.1 представлены первые шесть 
собственных резонансных частот, причем пятая частота для тела без 
трещины соответствует поперечному растяжению-сжатию в собственной 
форме колебаний вдоль линии стыковки материалов. 

Таблица 1 


Без дефектов, кГц Краевая задача 11, кГц 


2,851 1,721 
9,194 7,946 
9,858 9,575 


19,01 17.51 
22,32 27,18 
ПИТ У 24,16 27,88 





Решение ГИУ (21) проводилось на основе идей метода граничных 
элементов. При его дискретизации использовалась кусочно-постоянная и 
непрерывная кусочно-линейная аппроксимация неизвестных. Численные 
эксперименты показали, что более предпочтительной является первая 
схема. Решение дискретного аналога уравнения (21) проводилось методом 
регуляризации А.Н.Тихонова [7], параметр регуляризации в случае п= 40 


был равен 2 107“ (17 - число точек, в которых проводились "измерения"). 
На рис.3 изображены напряжения б ..(х,0) 10`° н/м? (кривые 1) и 


0: (х,0) 10° н/м? (кривые 2), найденные из решения обратной задачи, 


причем на рис.3,а прямая задача не содержала трещин, а рис.3,6 - 
соответствует одной трещине К. В численном эксперименте изучали 
влияние количества И входной информации на точность реконструкции 
трещин. 
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На рис.4 изображены напряжения с ..(х,0) 10`° н/м? (кривые 1, 2) 


иб..(х,0) 10° н/м? (кривые 3, 4), причем кривые 2 и 4 представляют 
конечно-элементное решение прямой задачи с двумя трещинами КЕ и ММ, а 
кривые 1, 3 найдены из решения обратной задачи при и-= 40. Как видно из 
рис.3 и 4, концы трещин легко идентифицируются по экстремальным 
значениям напряжений при решении обратной задачи. 
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В табл.2 представлена относительная погрешность (в процентах) 
нахождения координат концов трещин хх, х,, хи, ху, когда шаг между 


точками "измерений" составлял 0,0025 м и они располагались в центре 
стороны ВС. 


Таблица 2 





Во втором примере идентифицируется интерфейсная трещина в 
свободной биморфной (Ре, Си) пластине (0.3х0.04 м). В табл.3 
представлены первые три собственных частоты и формы колебаний. Из 
этой таблицы следует, что частота в окрестности первой изгибной формы 
не подходит для эффективной реконструкции, при использовании второй 
собственной частоты главную роль играет горизонтальное смещение 
берегов трещины, а на третьей собственной частоте раскрытие трещины 
происходит по вертикали. 

Таблица 3 


Собственная частота в кГц Собственная форма колебаний 
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1.82 Распределение 
вертикального 
смещения 





2.57 Распределение 
горизонтального 
смещения 








5.97 Распределение 
вертикального 
смещения 











Для реконструкции трещины использовалось уравнение (17) при 
его дискретизации неизвестных функций - скачков перемещений берегов 
трещины с помощью 10 кусочно-постоянных элементов. Так, на рис.5 
представлены найденные скачки горизонтального смещения (слева) в 
окрестности второй собственной частоты (16000 рад/с) и вертикального 
смещения (справа) в окрестности третьей собственной частоты. 





























10 10 


Рис.5 


Проведенные численные эксперименты показали, что точность 
определения координат концов трещин в значительной мере зависит от 
расположения участка "измерения". 

Авторы благодарят сотрудников кафедры математического 
моделирования ЮФУ - коллег по разработке пакета АСЕГАМ, результаты 
которых использовались при проведении расчетов. 


Библиографический список 

1. Ваппоуг Т, Веп АБЧа А., Заоца М. А зетгЕ-ехрИс а!догиИт 
Гог Че гесопыгисНоп оЁ ЗО рапаг сгасК$. // Тпуегзе РгоШет. 1997. \.13 
—Р. 899-917. 


193 


Физико-математические науки 








2. Ватульян А.О. Определение ориентации плоских трещин в 
упругом теле / А.О. Ватульян, А.Н. Соловьев // Теорет. и прикладная 
механика. - 2003. — Вып. 37. - С. 141-145. 

3. Соловьев А.Н. Идентификация интерфейсной трещины. 
Современные проблемы МСС: тр. \УПШ международ. науч. конф. 
[ А.Н. Соловьев. - Ростов н/Д. - 2002. - Т.1. — С 163-169. 

4. Ватульян А.О. Об одном классе граничных задач в динамической 
теории упругости / А.О. Ватульян, И.И. Ворович, А.Н. Соловьев // ПММ. 
— 2000. -—Т. 64, В. 3. - С. 373-380. 

5. Меж \\/., Апага Н., 5сппаск Е. Ап АКегпайпа Кегафуе адогНт Гог 
{Не гесоп&гисйоп оЁ ищегпа! сгасК$ т а {Игее-дитепз!опа! зойа Боду. // 
Гпуегзе РгоШет. 2001. \/.17. Р. 1957-1975. 

6. Козлов В.А. Итерационный метод решения задачи Коши для 
эллиптических уравнений / В.А. Козлов, В.Г. Мазья, А.В. Фомин // ВМ и МФ. 
— 1991.-Т. 31. - С. 45-52. 

7. Тихонов А.Н. Методы решения некорректных задач 
А.Н. Тихонов, В.Я. Арсенин. -— М.: Наука, 1979. - 285 с. 

8. Ватульян А.О. Некоторые полуявные алгоритмы реконструкции 
интерфейсных трещин / А.О. Ватульян, А.Н. Соловьев // Изв. вузов. 
Сев.-Кавк. регион. Естеств. науки. - 2003. — № 3. - С. 20-24. 

9. Новацкий В. Теория упругости / В. Новацкий. - М.: Мир, 1975. 
— 872 с. 

10. Белоконь А.В. Блочные схемы метода конечных элементов для 
динамических задач акустоэлектроупругости / А.В. Белоконь, В.А. Еремеев, 
А.В. Наседкин, А.Н. Соловьев // ПММ. - 2000. —Т. 64. — Вып. 3. - С.381-393. 

11. Наседкин А.В. Разработка визуализации и управления в 
конечно-элементном комплексе АСЕТАМ / А.В. Белоконь, А.Л. Петушков, 
А.Н. Соловьев // Зимняя школа по механике сплошных сред (тринадцатая). 
Школа молодых ученых по механике сплошных сред: тез. докл. 
— Екатеринбург: УрО РАН, 2003. - С. 274. 


Материал поступил в редакцию 11.02.09. 


А.М. ЗОГО\ТЕМ, А.5. 5РОСАКТМ, А.О. УАТОЕУАМ, $.М. 5НЕУТ$О\ 
ЕЕАМ/$ КЕСОМ$ТКОСТТОМ ТМ МУЕТТ-ГАУЕВ$ СОМРО$ТТЕ$ 


Треге аге оНеге4 тео4$ оЁ Таи $ гесоп$гисНоп тм $аНЯе сотрозКе$ - 
$(гаНЯсаНоп$ (е ме Масе ам). Тпуегзе деотеса| ргоБШетз$ о 
гесоп${гисНоп ог Яам/5 аге гедисе4 то а 5о!иНоп оЁ зедиепсе оЁ зу$етв оЁ {пе 
ВТЕ оатей оп {Ие Ба$5 оЁ гесргосйу Пеогет оЁ орегаЧоп$ аЁ е$а Неа 


194 


Вестник ДГТУ, 2009. Т9. №2(41) 








о$сШаЧоп$ оЁ ехрей затр/е$. 5оиЧоп ВТЕ мл! Бе сагйеа оиё оп {Пе Баз5 оРа 
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